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Abstract 

Some wireless channel models based on probability density 
functions have been proposed, including models based on 
Markov chains. However, they are limited to a certain 
number of perturbations. This paper proposes a 
methodology based in a Ku-band signal behavior which is 
classified in three cases: ascendant, descendent and constant. 
Samples are classified in these cases so their conditional 
probabilities are analyzed in order to find a fit to a 
probability density function and to extract its statistical 
parameters which are the model itself. From this information, 
a new signal was generated and their second order statistics 
are compared with the ones from the original signal to 
validate the created model. The analyzed signal was 
extracted from a measurement campaign done in Mexico. 
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Introduction 

A deep knowledge is required in the design and 
optimization of wireless communication systems. 
Models based on probability density functions (pdf) 
are not enough to characterize a dynamic channel 
because they only describe one type of phenomenon so 
other alternatives are analyzed. Many publications 
describe models created from measured samples 
where the signal is analyzed as a random variable and 
only one perturbation is characterized or it can only be 
applied to certain frequency or system. Other models 
are based in Markov chains which implies the 
combination of many phenomena however they are 
limited to a pre-established number of states or events 
during the transmission [(Barts, 1988), (Fontan, 2001), 
(Abouraddy, 2000), (Kattenbach, 2002), (Eberlein, 
2007)]. 

Another alternative for channel modeling is to analyze 
the variations in the level of the signal due to all 
phenomena (like rain, scintillation, etc.) which 



produces its fading. These variations are denominated 
as the dynamic of the wireless signal. The classification 
according to this dynamic is suggested by U.C. Fiebig 
in (Fiebig, U.-C, 1999) where samples which contain the 
attenuation of a signal are classified in three cases: 
ascending, descending and constant. An example of 
this classification is shown in FIG. 1. 

A wireless channel model for a set of samples from a 
Ku band signal is proposed in this paper. The model is 
based on works of (Fiebig, 1999) where a classification 
of the samples is done according to its dynamic, next a 
goodness of fit test is done to know which pdf fit best. 
The objective is to extract the statistical parameters of 
the fitted pdf for every value of level of the signal and 
case. These data are the model itself and a new time 
series can be generated from them. The validation of 
that process is done by comparing the second order 
statistics (average fading duration -AFD - and level 
crossing rate - LCR) of the original signal and the new 
generated signal. The document is organized as 
follows: the measurement campaign organized to 
extract the samples of the signal is described in section 
II; in section III the methodology is explained and it 
will be applied to data in order to obtain a new signal; 
results are explained in section IV; conclusions are 
given in section V. 
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FIG. 1 CLASSIFICATION OF SAMPLES ACCORDING TO ITS 
DYNAMIC 
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Measurement Campaign 

The signal whose behavior was studied comes from 
Intelsat 9 satellite which transmits the SKY Mexico's 
direct- to-home TV service (LyngSat, 2012). The central 
frequency of the analyzed signal is 12 GHz. The 
measurement campaign was done from August 25th to 
November 30th, 2011 at Tecnologico de Monterrey 
Campus Estado de Mexico. FIG. 2 shows the diagram 
of the measurement system used. 




FIG. 2 MEASUREMENT SYSTEM FOR THE MEASUREMENT 
CAMPAIGN 

Since rain is the principal source of attenuation in Ku 
band, quantity of rain was used as a base to organize 
the collected data. This parameter was taken from 
SEMARNAT data base whose sampling time is 10 min 
(Secretaria, 2010). According to this information, two 
periods were identified, one with a higher average of 
quantity of rain than the other. In this paper, all 
samples that corresponds to the higher quantity of rain 
were analyzed which covers from August 25 th to 
October 15 th , 2011 with average quantity of rain of 
0.0419 mm. The measured signal for this period is 
shown in FIG. 3 and FIG. 4 shows the samples of 
quantity of rain; for both FIG.s, min corresponds to 
the first taken sample for the analyzed period. 

In this type of data processing it is recommended to 
take only that samples that represents all the set of data, 
reducing in this way the quantity of information and its 
processing time (Husson, L., 2002). This process was 
done by using principal component analysis (PCA) 
which is a non-parametric method where the synthesis 
of the information is done for the purpose of reveling 
structures based on the variance and covariance of the 
analyzed data. The parameter to determine the most 
important information was based on Velicer algorithm 
and by using (1). 



W) 2 



(1) 



Where, 

f k , Sum of the square of the partial correlations in k. 



r£ , Correlation matrix R k . 

p„ Number of dimensions of the original set of 
data. 

Methodology 

The input data to begin the process are the results 
obtained after PCA where data are reduced and thus 
the time of processing too. Then, from the observation 
of the signal, three cases are identified: increasing or 
decreasing of it level or maybe stays near the same 
point. The cases are characterized as it is shown in (2). 
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FIG. 4 QUANTITY OF RAIN 

r{t) - r(t - 1) < Descendent 

r(t) -r(t-l)>0 Ascendant (2) 

r(t) - r(t - 1) = Constant 



Classification of the signal according to 
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its dynamic. 

r(t), 

r(t-l) 



by following (8). 



Current sample. 
Previous sample. 



Data are separated in groups according to (2) and each 
group is analyzed in pairs of samples in order to 
obtain the conditional probabilities of two consecutive 
events including the case where the same event 
appears. The equation for this step is shown in (3). 

Px\Y ^dynamic (0 ' dynamic (? ~ I)) = 

P{* = »W (0 1 Y = -D} = 

P{x = r dymmic (t ), Y = r dymmic (t - 1)] 



P^J dynamic (0» T dynamic ~ I)) 



(3) 



Where, 



dynamic > 



dynamic 



Current sample, 
(f - 1) , Previous sample. 



This methodology is based on Fiebig's work [6]. Also 
this process is similar when a Markov chain analysis is 
done (Diana Alejandra Sanchez-Salas, 2008) but in this 
case there are not a limited number of states. 

A Chi-square goodness of fit test with a confidence 
interval of 95% is done for each level value and case 
(ascending, descending and constant). The null 
hypotheses for each test are in (4) - (7). 



First null 
hypothesis: 

Second null 
hypothesis: 

Fourth null 
hypothesis: 

Fifth null 
hypothesis: 



H : data fits a Normal pdf (4) 



H : data fits a Weibull pdf (5) 



H : data fits a Rice pdf (6) 



H : data fits a Rayleigh pdf (7) 



The choice of these pdfs is because they are the most 
common used (Proakis, 2008). Once it is determined 
the pdf with best fit, three tables containing the 
statistical values with the pdf, one for each case 
(ascendant, descendent and constant) and for each 
possible level value of the signal in dB. Data from these 
tables are the input parameters for the time series 
generator. The generated samples are analyzed in pair 



^Isamples U t-\ U t 



(8) 



mples , Difference between two samples. 



Where, 
A 

nt-1, First or precedent sample. 

nt Second or current sample. 

Computed results from (8) are classified into three 
cases as it is shown in (9). 



dynamic 



Where, 



2 samples 



- 1 Constant 



Ao_, c > 1 Descendent 

1 samples 

A 9 , <1 Ascendant 

2 samples 



(9) 



dynamic 



, Classification of the generated signal according 
to its dynamic. 

^2 samples ' Difference between two samples. 

These three possibilities are the same as original cases 
of (2) so it is deduced this is an inverse process. Once 
the case was identified, a new sample is generated 
based on the statistical parameter that corresponds to 
the case and the current level of the signal. The final 
step is to validate the generated signal is statistically 
similar to the original signal through second order 
statistics (Perez Fontan, 2008). 

Results 

Ku-band Channel Model 

From PCA it was obtained that all samples 
corresponding to the first eight days contains the most 
important information and the firsts four days for 
condition 2. These are the samples that will be 
analyzed. The results after this classification are shown 
in FIG. 5. 

Next, samples were analyzed in pairs to compute the 
conditional probabilities of each event. Then, Chi 
square goodness of fit test for each case and level value 
in dB of the signal with confidence interval of 95% 
were applied. The null hypotheses are listed below. 

1. H 0a : Data fit a Normal pdf. 

2. H : Data fit a Weibull pdf. 
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3. H : Data fit a Rice pdf. 

4. H : Data fit a Rayleigh pdf. 

Results of those tests are shown in Tables 1 -3. Since 
the conditions of this type of test is there must be at 
least 5 elements in order to obtain reliable parameters, 
this test was applied to some level values of the signal 
because there were not enough information. This is the 
reason why not all possible values are shown. 

Tables 1-3 shows a better fit to a Weibull and Rice pdf, 
for the other tests there are not enough information to 
confirm the fit. In order to characterize the signal, 
Weibull pdf was chosen as the best fit. The parameters 
for each level of the signal and case are shown in 
Tables 4-6. 

TABLE 1 RESULTS OF GOODNESS OF FIT TEST FOR 
DESCENDENT CASE 
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+ Constant 
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[dB] 
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FIG. 5 CLASSIFICATION OF THE DYNAMIC OF THE MEASURED 
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TABLE 2 RESULTS OF GOODNESS OF FIT TEST FOR 
ASCENDANT CASE 



Level 
[dB] 
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TABLE 3 RESULTS OF GOODNESS OF FIT TEST FOR CONSTANT 

CASE 



TABLE 4 STATISTICAL PARAMETERS FOR DESCENDENT CASE 
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TABLE 5 STATISTICAL PARAMETERS FOR ASCENDANT CASE TABLE 6 STATISTICAL PARAMETERS FOR CONSTANT CASE 



Level [dB] 


k 
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Generation of the Signal and Validation of the Model 

Validation of the created statistical model includes the 
generation of a new signal from extracted Weibull pdf 
parameters for each level in dB and case. The 
generated signal is shown in FIG. 6. A comparison 
through second order statistics of the original analyzed 
signal and the generated one was done. Results of that 
comparison are shown in FIG. 7 and 8. 




FIG. 6 GENERATED SIGNAL FROM EXTRACTED STATISTICAL 
PARAMETERS 
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FIG. 7 COMPARISON OF AFD BETWEEN ORIGINAL AND 
GENERATED SIGNAL 

Generated signal of FIG. 6 was created by using 
Weibull parameters of Tables 4-6 and at first sight it 
seems similar to the signal of FIG. 3. Second order 
statistics confirm this hypothesis since both curves are 
similar in shape and values, especially in level crossing 
rate of FIG. 8 where an on-off behavior around -38 dB 
is appreciated which is a characteristic of urban 



scenarios (Scalise, 2005) as the one where measures 
were taken; it also confirms that the values of the 
measured signal concentrates in two sets: from -30 to 
-37 dB and from -39 to -45 dB. For the average fading 
duration there is a higher point in -38 dB which 
indicated signal remains more frequently under that 
level and descends for the next levels but it preserves 
the ascendant tendency. The values in average fading 
curves of the original and generated signal are not 
exactly the same but they have the same shape and 
characteristics of a typical AFD curve (Perez Fontan, 
2008). In this way it is concluded that a similar 
statistical signal was created from the proposed 
channel model. 
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FIG. 8 COMPARISON OF LCR BETWEEN ORIGINAL AND 
GENERATED SIGNAL 

Conclusions 

A Ku-band signal was analyzed in this paper in order 
to create a wireless channel model. That signal was 
obtained from a measurement campaign done in 
Mexico which contains a TV signal. During the 
measurement campaign it was detected two different 
periods based on the quantity of rain. Since rain is the 
principal source of attenuation, samples corresponding 
to that period were used; the average quantity of rain 
was 0.0419 mm. The model was based on the analysis 
proposed by U.C. Fiebig in [6] where three types of 
movements were identified in the signal: ascendant, 
descendent and constant. From this classification, data 
were examined according to the cases and the amount 
of events by analyzing the current and precedent 
samples. After data was analyzed, the probability of 
each event is computed, and then a goodness of fit test 
is applied to determine if data (according to the case 
and attenuation value) fits a specific probability 
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density function so its parameters were computed. The 
parameters of the pdf according to the case and level 
values are the channel model itself because with these 
data a new time series can be generated. For the 
analyzed signal, data fit a Weibull pdf. The last phase 
of this model is validation of it which consists in the 
comparison of second order statistics: LCR and AFD 
were curves of original and generated signal were very 
similar concluding in this way those Weibull 
parameters for each case characterizes the wireless 
channel. The methodology presented in this paper 
could be used to obtain a wireless channel model of 
any band to analyze the behavior of the signal. The 
generation of a new signal is useful to simulate the 
behavior of the wireless channel in the design of new 
telecommunication systems or to analyze new fading 
mitigation techniques. 
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